Comparison of post-Newtonian templates for extreme mass ratio inspirals 
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Extreme mass ratio inspirals (EMRIs) , the inspirals of compact objects into supermassive black 
holes, are important gravitational wave sources for the Laser Interferometer Space Antenna (LISA). 
We study the performance of various post-Newtonian (PN) template families relative to the high 
precision numerical waveforms in the context of EMRI parameter estimation with LISA. Expressions 
for the time domain waveforms TaylorTl, TaylorT2, TaylorTS, Taylor T4 and TaylorEt are derived 
up to 22PN order, i.e 0{v'^'^) {v is the characteristic velocity of the binary) beyond the Newtonian 
term, for a test particle in a circular orbit around a Schwarzschild black hole. The phase difference 
between the above 22PN waveform families and numerical waveforms are evaluated during two-year 
inspirals for two prototypical EMRI systems with mass ratios 10~* and 10~^. We find that the 
dephases (in radians) for TaylorTl and TaylorT2, respectively, are about 10~^ (10~^) and 10~® 
(10^^) for mass ratio 10"* (10"^). This suggests that using 22PN TaylorTl or TaylorT2 waveforms 
for parameter estimation of EMRIs will result in accuracies comparable to numerical waveform 
accuracy for most of the LISA parameter space. On the other hand, from the dephase results, we 
find that TaylorT3, Taylor T4 and TaylorEt fare relatively poorly as one approaches the last stable 
orbit. This implies that, as for comparable mass binaries using the 3.5PN phase of waveforms, 
the 22PN TaylorT3 and TaylorEt approximants do not perform well enough for the EMRIs. The 
reason underlying the poor performance of TaylorT3, TaylorT4 and TaylorEt relative to TaylorTl 
and TaylorT2 is finally examined. 

PACS numbers: 04.25. Nx 



> 

i> 

in 

^' 

o 
en 



X 



I. INTRODUCTION 

The inspiral of a stellar mass compact object into a su- 
permassive black hole (SMBH) is one of the most promis- 
ing gravitational wave (GW) sources for space based de- 
tectors such as cLISAQ The compact object typically 
has a mass of the order of a few solar masses while 
the SMBHs detectable by eLISA arc in the mass range 
IO^Mq — IO'^Mq. As the mass ratio for these binaries 
is typically around 10^, these systems are called extreme 
mass ratio inspirals (EMRIs) . Gravitational waves from 
EMRIs can provide information about the parameters of 
the central black hole such as its spin, mass and details 
of its stellar surroundings while also facilitating strong 
field tests of general relativity (GR) [ll-Q- 

However, the gravitational wave signal is buried in a 
background of noise and the signal needs to be extracted 
using data analysis techniques such as matched filtering. 
Matched filtering requires accurate templates of the grav- 
itational waveform. For eLISA EMRI parameter estima- 
tion, the GW phase errors of the template with respect 
to the true signal should be less than 10 milliradians [j]. 
Considering eLISA is expected to detect 10-1000 EMRIs 
during its mission [ll-Q, the search for accurate wave- 
forms for these systems is justified. 



^ evolved Laser Interferometer Space Antenna, also known as NGO 
(New Gravitational- Wave Observatory). 



Post-Newtonian (PN) theory provides a method to pre- 
dict the gravitational waveform for the early phase of 
inspiraling compact binaries. However, since the PN 
approximation breaks down near the last stable orbit 
(LSO), numerical relativity (NR) waveforms are required 
beyond this point. Nevertheless, PN waveforms in the 
early inspiral phase are required to calibrate the NR 
waveforms because the computational cost for NR is very 
high. PN waveforms can be matched with the NR wave- 
forms in the late inspiral and the subsequent merger and 
ringdown phases [8|, |9| to provide a cheaper alternative 
to using NR for the complete inspiral. 

PN waveforms for non-spinning comparable mass bina- 
ries in quasi-circular orbits are known up to 3.5PN order 
beyond the Newtonian term |10l - [l3 |. Within the post- 
Newtonian formalism, several non-equivalent template 
families such as TaylorTl, TaylorT2, TaylorTS, Tay- 
lorT4, TaylorEt and TaylorF2, among others, are possi- 
ble. These 3.5PN template families were discussed exten- 
sively in ReL [l^. It is found that for M < 12 Mq, where 
M is the total mass of the binary, these 3.5PN tem- 
plate families except for TaylorTS and TaylorEt are 
equally good for the detection of gravitational waves us- 
ing ground based detectors. 

Recently, 22PN waveforms have been calculated for a 
test particle in a circular orbit around a Schwarzschild 
black hole [l3- It is shown that the 22PN gravitational 
waveforms achieve data analysis accuracies comparable 
to those resulting from high precision numerical wave- 



forms. In this paper we extend this study by calculat- 
ing the different template families mentioned above up 
to 22PN order using the 22PN energy flux derived in 
Ref. [13 . We then investigate the performance of these 
Taylor approximants by evaluating the phase difference 
between these approximants and the high precision nu- 
merical waveforms in Refs. |14l . |l5| over a two-year in- 
spiral for two systems, one in the early inspiral phase 
and the other in the late inspiral phase of the eLISA fre- 
quency band. We find that TaylorTl (which was also 
investigated in Ref. [1^) and TaylorT2 provide the best 
matches to numerical waveforms while the phase differ- 
ence increases by a few orders of magnitude for TaylorTS, 
TaylorT4 and Taylor Et. These investigations extend the 
results for comparable mass binaries in Ref. [l3| that Tay- 
lorTS and TaylorEt approximants are considerably dif- 
ferent from the others and perform relatively poorly. We 
also discuss why the performance of TaylorTS, TaylorT4 
and TaylorEt in the test particle limit becomes worse 
than the others. This may provide insights which should 
be kept in mind when one constructs new PN template 
families. 

This paper is organized as follows. In Sec. II we dis- 
cuss the various template families along with the relevant 
initial conditions and calculate these approximants upto 
22PN order. In Sec. Ill we evaluate the dephase between 
these different PN waveform approximants and high pre- 
cision numerical waveforms during two-year inspirals. In 
Sec. IV we summarize our main conclusions. Since the 
22PN Taylor approximants are too large to be shown in 
this paper, we only show them up to 4.5PN order. The 
22PN expressions for the approximants will be publicly 
available online. Throughout this paper we use units 
c = G= 1. 



II. THE POST-NEWTONIAN APPROXIMANTS 

Post-Newtonian approximation treats the early stages 
of adiabatic inspiral of compact binaries as a pcrturbativc 
model and expresses a binding energy, E(v) and a flux, 
J^(v) associated with the gravitational wave as a power 
series in v, where v = (ttMF)^/'^ is the characteristic 
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velocity, M is the total mass and F is the gravitational 
wave frequency of the binary. Here, adiabatic inspiral 
implies that the inspiral time scale is much larger than 
the orbital time scale. For restricted waveformsp under 
the adiabatic approximation, the standard energy bal- 
ance equation, dEtot/dt ~ —J- gives us the following pair 
of coupled differential equations for the orbital phasing 
formula Jiii]J|: 



dt 



V 

JI 



0, 



Tiv) 



dv 

'dt ' ME'{v) 



= 0. 



(2.1a) 



(2.1b) 



Here, E'{v) is the derivative of the binding energy 
with respect to the characteristic velocity, v. The bind- 
ing energy, E is related to the total energy, Etot by 
Etot=M{l + E). 

The phasing formula can also be expressed in the fol- 
lowing equivalent parametric form, where frof and 0rof 
are integration constants and Wref is an arbitrary refer- 
ence velocity. 



t{v) = t,^f + M 



(v) = (/)rcf + 



■-■f E'{v) 



dv, 



tr—rH-dv. 
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(2.2a) 



(2.2b) 



Recently, the 22PN order energy flux for EMRIs has 
been calculated [l^- For the extreme mass ratio bina- 
ries, E(v) is known exactly and can be expanded and 
truncated to any required PN order. We present for the 
convenience of the reader, expressions for 4PN E{v) and 
4.5PN J^{v) that are inputs needed to derive the 4.5PN 
results displayed explicitly in later sections for brevity of 
presentation. In these expressions, mi and 7712 are the 
masses of the test particle and the SMBH respectively, 
M — TOi -|- 1712 is the total mass, v ~ toiTO2/M^ is the 
symmetric mass ratio and 7 = 0.577216... is the Euler 
constant. 
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^ Restricted waveforms are obtained by retaining only the leading 
harmonic in the GW signal. For these waveforms the gravita- 



tional wave phase is twice the orbital phase. 
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Approximate waveforms are obtained by inserting the 
expressions for E{v) and T{v) at consistent PN order 
in the phasing formula - these waveforms are referred 
to as Taylor approximants. There are several ways of 
inserting these expressions in the phasing formula leading 
to different approximants such as TaylorTl, TaylorT4, 
TaylorT2, TaylorT3, TaylorEt and TaylorF2 ^. We 
have calculated these approximants up to 22PN order for 
the EMRI case. We shall now discuss these approximants 
while presenting our results up to 4.5PN. Method for 
calculating the phase of the gravitational waveform is 
left for the next section. The complete 22PN expressions 
will be available online. 
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In the above equations v = v^"^^' , but for simplicity we 
omit the superscript. The expressions for E{v) and T{v) 
arc to be truncated to consistent PN order to obtain the 
approximant of that order. This is followed for all the 
approximants in this section. 



A. TaylorTl 

The TaylorTl approximant is obtained by using the 
expressions for binding energy, E{v) and flux, J^{v) as 
they appear in Eq. (|2.3p and Eq. (|2.4p in the phasing 
formula, Eq. (|2.ip and solving the resulting equations in- 
volving the rational polynomial T{v)/ E' [v) numerically. 
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B. TaylorT4 



TaylorT4 goes one step further than TaylorTl by ex- 
panding the rational polynomial J-{v)/E'{v) and trun- 
cating it to the required PN order. The characteristic 
velocity, v^'^'^\t) = v{t) at 4.5PN is given for TaylorT4 
by, 
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Similar to TaylorTl, Eq. (|2.1ap gives the evolution of 
orbital phase for TaylorT4. 



truncated to required PN order. Upon integration we ob- 
tain the following equations for (/)(w) and t{v) at 4.5PN 
order: 



C. TaylorT2 

TaylorT2 uses the parametric form of the phasing for- 
mula, Eq. ()2.2|) . The ratio E'{v)/F{v) is expanded and 
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Here trcf and 0rcf arc integration constants, iref is fixed 
by setting t = when v ~ vq, the initial velocity. 



D. TaylorT3 

To get the TaylorT3 approximant, the expression for 
t{v) generated in TaylorT2 is inverted to get v{t). This 
is then used to obtain (j){t) = (j){v{t)). TaylorTS also 
gives the instantaneous gravitational wave frequency, F 
hy F = dtfi/indt) = v^/^ttM). TaylorT3 approximant at 
4.5PN order is given by, 
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where e(t) is given by 6* = [iy{t,cf - t)/{5M)]-^^^. Given Fq, by Fq = v^/{-kM). To find tref, one solves Eq. ^^Ih^ 
an initial velocity vq^ one can find the initial frequency 



at i = and i^ = Fn 



E. TaylorEt 



From the phasing formula Eq. (|2.1a[) and Eq. p.9|) we 
get an expression for the evolution of phase in terms of ^ 
(cf. Eq. (j2.11a|) ). Under the new variable, C, Eq. (|2.1bp 
of the phasing formula transforms to 



TaylorEt is expressed as a power series of a new func- 
tion, ^ = —2E/v. Eq. (|2.3|) for E{v) can be expressed 
in terms oi x = v^ to get (^{x). This is then inverted to 
obtain x(C): 
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TaylorEt approximant is, essentially, the gravitational 
wave phasing equations expressed in terms of C. At 
4.5PN order, it is given by. 
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As in the case of TaylorTS, we can find Fq, given vq. 
Noting that F = d(j)/{Trdt), initial conditions for TaylorEt 
can be setup by solving Eq. (|2.11a[) for Co by setting the 
left hand side to ttFq. 



I 

tional wave frequency, F(t) becomes equal to the Fourier 
variable, /. In the adiabatic approximation (where 
Eqs. (|2.2ap and (|2.2bp hold) the values of i/ and ■(/'/(*/) 
arc given by. 



F. TaylorF2 
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TaylorF2 is a Fourier domain approximant based on 
the stationary phase approximation (SPA). Under the 
SPA. the waveform in the Fourier domain is expressed as 
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where tf is the saddle point, defined by solving for t 
when dtpf{t)/dt ~ 0, i.e the time, tf when the gravita- 
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(2.13b) 
where Vj = {irMfy/^. 

Using expressions of energy and fiux and expanding 
the ratio E'{v)/F{v) in Eq. (j2.13p to consistent PN order 
leads to an expression which can be integrated explicitly 
resulting in the TaylorF2 approximant. The phase of the 
Fourier domain waveform up to 4.5PN order is given by. 
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where tc and 0c can be chosen arbitrarily and v ~ 
(7rM/)i/3. 

The behavior of TaylorF2 has already been investi- 
gated up to 3.5PN order for comparable mass binaries 
[l0|-[l2| . One must keep in mind that the stationary phase 
approximation, on which TaylorF2 is based, is valid only 
up to 4.5PN order [H, [T^. Thus, beyond 4.5PN the 
Fourier transform of the waveform has correction terms 
for the stationary phase approximation to the Fourier 
transform in Eq. (|2.12p . However, as a start, in this pa- 
per, we have obtained the TaylorF2 approximant upto 
22PN order by assuming Eq. (|2.12p is valid even be- 
yond 4.5PN. Further studies, computing terms beyond 
the leading jl(]| - [l2| are needed to look for a good fre- 
quency domain approximant at higher PN orders and 
will be investigated in the future. 



III. COMPARISON WITH NUMERICAL 
RESULTS 



To investigate the behavior of different analytical PN 
families described in Sec. |lT] we calculate the phase of 
the gravitational wave signal during a two-year quasi- 
circular inspiral of two systems of binaries called System-I 
and System-II as considered in [IJ, [l8l - l20 |. We com- 
pare this phase with the phase calculated using high 
precision numerical waveforms, the difference between 
these phases is called the dephase. System-I has masses 
(mi, 7712) = (1O,1O^)M0 with mi/m2 = 10"'', it in- 
spirals from rjn — 29Af to rfjn ~ 16M during a two- 
year period, with gravitational wave frequencies in the 
range /gw S [4 x 10"'^, 10~^]Hz. System-II has masses 
(toi, 771,2) = (10, 1O^)M0 with 7711/7712 = 10~^, it inspi- 
rals from rin ~ IIM to ran — 6Af (LSO) during a two- 
year period, with gravitational wave frequencies in the 
range /gw G [1.8 x 10"^, 4.4 x 10~^]IIz. In the frequency 
band of eLISA, System-I corresponds to the early inspi- 
ral phase of an EMRI while System-II corresponds to the 
late inspiral phase. Note that the phases shown in the 
figures of this section are gravitational wave phases which 
are twice the orbital phases. 

Matched filtering can give signal to noise ratios (SNRs) 
of up to p ^ 100 [m, 0] for the strongest EMRI sig- 
nals detectable by eLISA. This means eLISA can detect 
bases up to an accuracy of order 1/p ~ 10 milliradians 
Therefore, while considering the dephase results of 
this paper, we expect PN waveforms to have accuracies 
comparable to those provided by numerical waveforms if 
the dephase is less than 10~^ radians. 

The numerical waveforms we use are based on those in 
Refs. [IJ, [lal- The truncation of the Z-mode limits the 
accuracy of the numerical calculations. We use the same 
data generated for Ref. [lj|, which is based on I = 25 
calculations and gives relative error better than 10^^^ up 
to the LSO. 



Dephase between TaylorTl and numerical 
results 



The dephase between TaylorTl and numerical wave- 
forms (cf. Fig. [1]) was shown in Ref. [13]. We present 
the results here for comparison^ For System-I (II), the 
absolute values of the dephasing between the TaylorTl 
waveforms and the numerical waveforms after the two- 
year inspiral are about 7 x 10^(8 x 10^), 7 x 10^^(8), 7 x 
10-6(5 X 10"i),8 X 10-9(3 X 10-2) and 10-^(5 x 10-^) 
radians for 5.5PN, lOPN, 14PN, 18PN and 22PN respec- 
tively. It is also suggested in Ref. [IJ] that using 22PN 
TaylorTl waveforms for EMRIs will result in accuracy 
of data analysis comparable to those resulting from high 
precision numerical waveforms as the dephase is less than 
10-2 radians for most of the parameter space of eLISA. 
We also note that lOPN TaylorTl waveforms may be 
comparable in accuracy of data analysis to numerical 
waveforms for System-I. 

We now extend this study by investigating the behav- 
ior of other PN Taylor families by evaluating their de- 
phases during the same inspiral. 



Dephase between TaylorT4 and numerical 

results 



Calculation of phase of TaylorT4 is very similar to Tay- 
lorTl and numerical waveforms. We use the relation 



m 



{d<j)/dt') dt' = 



„ {dv'/dV) 



dv'. (3.1) 



Here, dcjy/dt is given by Eq. (|2.1ap . vq is the velocity at the 
starting of the two-year inspiral and v{t) can be obtained 

by solving the equation t — J l/{dv/dt) dv = for a 
given time t. For TaylorT4 dv/dt is given by Eq. (|2.6p 
while for numerical waveforms dy/dt is obtained from the 
solution of Teukolsky equation J13l - ll5l | . 
The dephase results are shown in Fig. 



2(a) and Fig. 



2(b) For System-I (II), the absolute values of the de- 
phasing between the TaylorT4 waveforms and the nu- 
merical waveforms after the two-year inspiral are about 
2 X 10^(3 X lO'^), 7 X 10-1(7 x 10^), 6 x 10-^(2 x 10^), 5 x 
10-5(8 X 10^) and 5 x 10-^(3 x 10^) radians for 5.5PN, 
lOPN, 14PN, 18PN and 22PN respectively. This suggests 
that 14PN or higher order TaylorT4 waveforms are com- 
parable in accuracy of data analysis to numerical wave- 
forms for the early inspiral phase (System-I) but the ac- 
curacy is low for the late inspiral phase (System-II), par- 
ticularly, near the LSO. 

By comparing with Fig. [U we see that the dephase for 
TaylorT4 is a few orders of magnitudes worse than that 



^ Fig. [T]is slightly different from the dephase results of Ref. [l3l |. 
This is because in Ref. Il3ll the TaylorTl phase was calculated 
without expanding dE/dv in v. 
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(b)Systcm-II for TaylorTl 



FIG. 1. Absolute value of the dephase between TaylorTl FN waveforms and numerical waveforms during two-year inspirals 
as a function of time in months. The left panel shows the dephase for System-I having masses (mi,m2) = (10, 1O^)M0 for 
Tin — 29M to rfln — 16M sweeping GW frequencies in the range /gw G [4 x 10""^, 10~^]Hz. The right panel shows the dephase 
for System-II having masses (mi,m2) = (10, 1O^)M0 for Hn — IIM to rfln ~ 6M (LSO) sweeping GW frequencies in the range 
/gw G [1.8 X 10~^,4.4 X 10~'']Hz. System-I (System-II) corresponds to the early (late) inspiral phase of the eLISA frequency 
band. Note that the dephase between 18FN (22PN) TaylorTl waveforms and numerical waveforms at the end of the two-year 
inspiral for System-I is about 8 x 10~^ (10~^) radians, which falls below the lowest value of dephase in the left panel. 



for TaylorTl. This can be explained as follows. In the 
test particle limit, dE / dv, given by 



dE 
dv 



(l-6t;') 
(l-3f 2)3/2' 



(3.2) 



goes to zero as one approaches the LSO at w = l/-\/6. 
Therefore the series expansion of (dE/dv)^^ converges 
very slowly around the LSO. Noting that TaylorTl and 
TaylorT4 differ in whether or not the series expansion 
of {dE/dv)~^ is performed in obtaining dv/dt, one can 
expect that TaylorTl is more accurate than TaylorT4 
near the LSO. This suggests that factorization to avoid 
a pole at the LSO leads to improvement in the accuracy 
of dv/dt. We note that in Refs. [23, [2J] factorization is 
performed for the energy flux, T{v) to deal with a pole 
at the light ring, v = l/v3. Thus, one may also have to 
factorize the pole at the light ring when considering the 
case beyond the LSO. 



Dephase bet^veen TaylorT2 and numerical 
results 



TaylorT2 expresses the orbital phase, ^(w) and the 
time, t{v) as functions of v. For a given time T, we 
solve the equation T — t{v) = to get the velocity, v{T). 
Given the velocity one can compute the phase as 



<j>it) = ^(vit)). 



irct in Eq. (|2.7ap is chosen such that (j){ t ~ ) = 0. 
The dephase results are shown in Fig. 2(c) and Fig. 



|2(d)[ For System-I (II), the absolute values of the de- 
phasing between the TaylorT2 waveforms and the nu- 
merical waveforms after the two-year inspiral arc about 
9 X lOi(lO^), 4x 10-2(6 x 102), 5 x 10-'^(5), 5 x 10-^(10-2) 



and 10-9(8x10-") radians for 5.5PN, lOPN, 14PN, 18PN 
and 22PN respectively. Therefore for TaylorT2, lOPN 
(18PN) waveforms may provide accuracies comparable 
to those provided by numerical waveforms for System-I 
(Systcm-II). 

As can be seen by comparing with Fig. [l] the dephase 
of TaylorT2 is comparable or lesser than TaylorTl during 
the inspirals. However, one needs to keep in mind that 
for calculating the phase from TaylorT2 approximant a 
pair of transcendental equations needs to be solved which 
can be very time consuming and expensive. 



D. Dephase between TaylorEt and numerical 

results 



TaylorEt expresses dcji/dt and dQ/dt as power series 
expansions of C = —lEjv. For a given time T, we solve 
the equation T - f^^^ ll{dCldt)dC == to get C{T), 



where Co and C,{T) arc the values of C at t 
T, respectively. 

The phase can now be evaluated as 



and time 



0(t) ^ / {d(j)/dt')dt' 



at) 



{d(p/dt') 



dC, 



(3.4) 



where {d(f)/dt') and [dC/dt') are given in Eq. ([2lT|) . 
The dephase results are shown in Fig. 2(e) and Fig. 



(3.3) HSI 



For System-I (II), the absolute values of the de- 
phasing between the TaylorEt waveforms and the nu- 
merical waveforms after the two-year inspiral are about 
103(6x 10'*),9(2x 104),3x 10-i(9x 103),8x 10-3(7x 10^) 
and 4 X 10-''(4 x 10^) radians for 5.5PN, lOPN, 14PN, 
18PN and 22PN respectively This suggests that 18PN 
and 22PN TaylorEt waveforms arc comparable in accu- 
racy of data analysis to numerical waveforms for the early 
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FIG. 2. Absolute value of the dephase between difTerent PN waveforms and numerical waveforms during two-year inspirals as a 
function of time in months. The left (right) panel shows the dephase for System-I (System-II). Note that the dephase between 
18PN (22PN) TaylorT2 waveforms and numerical waveforms at the end of the two-year inspiral for System-I is about 5 x 10"'' 
(10"") radians, which falls below the lowest value of dephase shown in Fig. 2(c) If the dephase is less than 10 milliradians, 



the PN waveforms will provide data analysis accuracies comparable to those provided by the numerical waveforms. Therefore 
we expect that 18PN waveforms for TaylorT2 are comparable in accuracy of data analysis to numerical waveforms for most of 
the EMRI parameter space of eLISA. 14PN and 18PN waveforms are required for TaylorT4 and TaylorEt respectively, to be 
comparable to numerical waveforms even in the early inspiral phase (System-I). The accuracy of TaylorT4 and TaylorEt for 
System-II (which goes up to the LSO) is not comparable to numerical waveform accuracy. 



inspiral phase (Systcni-I) of cLISA frequency band but 
the accuracy is low for the late inspiral phase (System-II). 
As in the case of TaylorT4 in Sec. lIIIBl the reason that 
the performance of TaylorEt is much worse than that of 
TaylorTl can be related to the poor convergence of the 
series expansion of {dE/dv)~^ or dv/dt around the LSO. 



Solving C = —2E/i' iterativelyO the new variable C in 



* One can solve f = —2Elv to obtain the explicit expression for 
i)(^) without performing a series expansion in terms of v. But 



TaylorEt can be related to v as in Eq. (|2.9|) . One can also 
derive the same relation using v{Q = J{d(^' /dv)^^ dC,' . 
Noting {dC/dv)~^ = -v [dE/dv)"^ /2, we see that the 
integrand of v{C} contains a pole at the LSO. Then, one 
may expect that v as a series expansion in terms of C, does 
not converge very well around the LSO. Hence functions 
of C in TaylorEt, which arc computed by using a scries 
expansion of v(C), will not converge well. 



E. Dephase between TaylorTS and numerical 
results 

TaylorTS gives the orbital phase, <j){t) and the instan- 
taneous gravitational wave frequency, F{t) as functions 
of 9{t), which is a function of time. For any given time 
t, we can find the phase as 

0(t) - 4>{e{t)). (3.5) 

As in TaylorT2, 0rof is chosen such that (f>{t = 0) = 0. 
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FIG. 3. Absolute value of the dephase between TaylorTS PN 
waveforms and numerical waveforms during two-year inspirals 
for System-I as a function of time in months. TaylorTS is 
found to behave poorly for System-II as it goes up to the last 
stable orbit (LSO). We expect that 22PN TaylorTS waveforms 
are required to get data analysis accuracies comparable to 
those provided by numerical waveforms for System-L 



The dephase results for System-I are shown in Fig. 
[31 For Systcm-I, the absolute values of the dephas- 
ing between the TaylorTS waveforms and the numer- 
ical waveforms after the two-year inspiral arc about 
2 X 10^ 9 X 10\ 7 X 10~\ 7 x 10"^ and 3 x 10"^ radians 
for 5.5PN, lOPN, 14PN, 18PN and 22PN respectively. 
Therefore we see that 22PN TaylorTS waveform is re- 
quired to get data analysis accuracies comparable to nu- 
merical waveform accuracy for System-I. 

TaylorTS approximant is not accurate in the case of 
System-II as it goes up to the LSO (v ~ l/\/6). Even for 
System-I we find the dephase is a few orders of magnitude 



the explicit expression for i'(^) may not be useful to implement 
since it contains the square root of polynomial functions of C,. 



higher than that for TaylorTl. One also finds that the 
value of F^'^^if) becomes very large for 9(t) > 0.67 and 
for System-II one cannot find tref consistent with the one 
derived for TaylorT2. The reason for this is again sim- 
ilar to the reason for the poor behavior of TaylorTd as 
compared to TaylorTl. In the TaylorTS approximation, 
one derives v{t) as a series of t, i.e. 9, by iteratively in- 
verting t{v) in TaylorT2. v{t) as a series of t can also be 
derived using v{t) = J{dv/dt')dt' . Since the integrand 
of v{t), dv/dt, has a pole at LSO, one will expect that 
v{t) as a series of t does not converge well around LSO. 
Thus, functions oi 9 in TaylorTS, computed using v{t) in 
a series of t, will not converge well. 



IV. CONCLUSIONS 

Using the 22PN expression for flux, J^{v), derived in 
Ref. [ii], we calculated the TaylorTl, TaylorT2, Tay- 
lorTS, TaylorT4 and TaylorEt approximants up to 22PN 
order for a test particle in a circular orbit around a 
Schwarzschild black hole. We evaluated the performance 
of the PN waveforms by calculating the gravitational 
wave phase predicted for two EMRI systems, System-I 
(7711/7712 = 10~^) and System-II {mi/m2 = 10~^) dur- 
ing two-year inspirals. System-I (System-II) corresponds 
to the early (late) inspiral phase of the eLISA frequency 
band. The phase predicted by PN waveforms is com- 
pared with the phase predicted by high precision numer- 
ical waveforms for the same inspirals. For accurate eLISA 
EMRI parameter estimation with these PN waveforms, 
we need the difference of the phases, the dephase, to be 
less than 10~^ radians [J]. 

We found that the dephase between the 22PN wave- 
forms and numerical waveforms after a two-year inspiral 
for System-II is smaller than 10~^ and 10""^ radians for 
TaylorTl and TaylorT2 respectively. Therefore we ex- 
pect that these 22PN waveforms can be used to attain 
data analysis accuracies comparable to those provided by 
high precision numerical waveforms for most of the pa- 
rameter space of EMRIs. Moreover, for the early inspiral 
phase, lOPN waveforms for TaylorTl and TaylorT2 may 
be used for data analysis. 

However, the dephase of TaylorT4, TaylorEt and Tay- 
lorTS waveforms goes to values higher than 10^ radi- 
ans for System-II. This suggests that these approximants 
cannot be used for data analysis of late inspirals. We note 
that our results reinforce investigations in Ref. [10] that 
TaylorEt and TaylorTS are recommended not to be used 
for data analysis of comparable mass binaries. 

For Systcm-I we found that 14PN or higher PN or- 
der waveforms are required for TaylorT4, TaylorEt and 
TaylorTS to achieve comparable results in data analysis 
to using high precision numerical waveforms even in the 
early inspiral phase. We also found that the reason the 
dephases of TaylorT4, TaylorEt and TaylorTS waveforms 
are much larger than that of TaylorTl and TaylorT2 may 
be related to the fact that {dE/dv)~^ has a pole at the 



10 



LSO. This suggests that when construeting templates for 
coaleseing compact binaries, approximants avoiding the 
pole at the LSO arising from the energy function by fac- 
torizing it as in TaylorTl or introducing a new variable 
which cancels the pole may perform better. We hope that 
these studies also provide insights to construct more ef- 
ficient templates for coalescing compact binaries in the 
comparable mass case. 
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